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High-order harmonic generation (HHG) is investigated for HJ and its isotopomers under seven- 
and ten-cycle trapezoidal laser pulses at 800 nm wavelength and 7=4x10^^ W/cm^ intensity. We 
solved numerically full-dimensional electronic time-dependent Schrodinger equation with and with¬ 
out the Born-Oppenheimer approximation. We show that the HHG at the falling edge of a trape¬ 
zoidal laser pulse can result in redshift and complexity on the total HHG spectrum which can be 
removed by considering different laser pulse duration and nuclear motion not possible for sin^ and 
Gaussian laser pulses. We resolve the redshifts and complexities of the HHG spectra into different 
electronic and vibrational states and their interferences. 

PACS numbers: 42.65.Ky, 42.65.Re, 42.50.Hz, 33.80.Rv 


I. INTRODUCTION 

The non-perturbative interaction of intense laser pulses 
with atoms and molecules gives rise to the generation 
of high energetic photons well known as high-order har¬ 
monic generation (HHG) [ 1 , 2 ]. The HHG is well under¬ 
stood by a three-step model proposed by Corkum [3] and 
extended by Lewenstein et al. [4]. First, an electron 
tunnels into continuum from suppressed potential cre¬ 
ated with combinations of system’s Coulomb potential 
and laser pulse field. Then, the electron oscillates in the 
laser field and moves away from the ion core, and after 
the field reverses, it comes back to the core. Finally, the 
electron may recombine with its parent ion, leading to 
photon emission of high frequencies which are multiples 
of the frequency of the driving laser. The HHG is used 
to produce single isolated or trains of attosecond laser 
pulses which are necessary for the real-time observation 
of electronic dynamics [5]. The emitted photons in the 
HHG possess also time-dependent structural information 
on sub-femtosecond time scales which can be used to ob¬ 
tain the details of a system under investigation [6-10]. For 
example, from amplitude [ 11 ] and frequency [ 12 ] modu¬ 
lations in the HHG spectrum, molecular internuclear dis¬ 
tance can be obtained. 

The HHG process in molecules is more complex than 
that in atoms due to nuclear motion, two-center interfer¬ 
ence and different orientations of a molecule with respect 
to the laser field. We mention some recent works on the 
effect of nuclear motion on the HHG as follows. Experi¬ 
ments on H 2 and D 2 [7,13-14], and CH 4 and CD 4 [7] show 
that the HHG is suppressed more for the lighter isotopes 
due to faster nuclear motion. The nuclear motion can 
shorten the length of attosecond laser pulse trains [15] and 
can produce an isolated attosecond laser pulse because of 
long-trajectory suppression [16,17]. It is reported that 
even-order harmonics are produced when superposition 
of two electronic states is formed for large internuclear 
distances as a result of the system’s symmetry breaking 
[18]. It is also shown that in the presence of nuclear mo¬ 
tion, like in HJ and H 2 , the energy difference between the 


two lowest bonding and anti-bonding electronic states de¬ 
creases for large internuclear distances, thus anti-bonding 
electronic state is populated and this electronic state also 
contributes to the HHG [19]. The two-center interfer¬ 
ence minimum disappears at large internuclear distances 
where recombination to the first excited (anti-bonding) 
state becomes important [20]. The suppression of long 
trajectories and the HHG yield at different intensities due 
to nuclear motion are also reported [ 21 ]. 

By considering nuclear movement and laser pulse du¬ 
ration, Bian et al. observed only the redshift of HHG 
spectrum for a sin^ laser pulse [12]. For laser pulses with 
sin^ and Gaussian envelopes having only rising and falling 
parts, the blueshift or redshift of harmonics relative to 
odd harmonic orders occur inevitably at the rising and 
falling parts, respectively (for more details, see ref. [ 12 ] 
and references therein). In practice, it is more difficult 
to observe only the redshift than the blueshift. At some 
conditions for sin^ and Gaussian laser pulses, depending 
on pulse duration, type of the isotope and laser pulse in¬ 
tensity, both blueshift and redshift of harmonics can be 
present, leading to broadening of the harmonics if the 
HHG happens equally at both rising and falling edges 
of the laser pulse [ 12 ]. Therefore, the HHG spectrum 
certainly changes in the case of sin^ and Gaussian laser 
pulses. 

In this work, we investigate the effect of nuclear move¬ 
ment on the HHG spectrum under trapezoidal laser 
pulses, mainly focused on the laser falling part. We ad¬ 
dress some questions for different isotopes under trape¬ 
zoidal laser pulses not already introduced and answered: 
Would the rising and falling parts of a trapezoidal laser 
pulse affect the HHG? To what extent the effects would 
be dependent on the type of the isotope? What molecu¬ 
lar electronic and vibrational components contribute into 
the effects? Finally, how could we control the shape of 
HHG spectrum due to the corresponding effects? 

For the trapezoidal laser pulses, we show that the red¬ 
shift of harmonics may occur depending on the isotope 
and pulse duration. It is shown that the effect of laser 
falling part on the HHG spectrum could be even re¬ 
moved and therefore controlled for a trapezoidal laser 
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pulse which is not possible for sin^ and Gaussian laser 
pulses as mentioned above. 

We resolve the redshifts and complexities of the HHG 
spectra into different electronic and vibrational states and 
their interferences. Similar decomposing the HHG into 
different electronic [19,20] and vibrational [18] states has 
been reported but here we apply them to resolve the ob¬ 
served redshifts. We show that not only individual elec¬ 
tronic and vibrational states are involved in the redshifts 
but also their interferences are significant and should be 
taken into account. To our knowledge, resolving the red¬ 
shifts into their components have not been reported, even 
for sin^ and Gaussian laser pulses. To see whether the 
nuclear motion affect the HHG, we performed the calcu¬ 
lations for both fixed and freed nuclei. For the latter, 
we used different isotopomers in order to have a better 
understanding of the effects of the nuclear movement on 
the observed complexities in the HHG. 

Here, full-dimensional electronic time-dependent 
Schrodinger equation (TDSE) beyond the Born- 
Oppenheimer approximation (NBO) is solved numeri¬ 
cally for H^, and (X is a virtual isotope of H 
being 10 times heavier). The full-dimensional electronic 
TDSE within the Born-Oppenheimer approximation 
(BO) for HJ, which is indicated throughout the article 
by HJ(BO), is also solved to compare with NBO results. 
The equilibrium internuclear distance within BO is set 
to Req = 1.96 a.u. Galculations have been done with 
seven- and ten-cycle laser pulses of 800 nm wavelength 
with I =4 xlO^^ W/cm^ intensity. We suppose that 
molecular ions are aligned in such a way that their 
internuclear-distance axis is parallel to the laser polar¬ 
ization direction. Nowadays the molecular alignment 
in strong laser fields is readily feasible and it is applied 
in several experiments [6,9,22-27]. We use atomic units 
throughout the article unless stated otherwise. 


II. COMPUTATIONAL METHODS 


The time-dependent Schrodinger equation for H J (D J) 
with electron cylindrical coordinate (z^p) with respect 
to the molecular center of mass and internuclear distance 
i?, for both 2 : and R parallel to the laser polarization 
direction, can be read (after elimination of the center-of- 
mass motion and ignoring molecular rotation) as [28-29] 


dt 


H{z,p, R;t)'tp{z,p,R;t). 


( 1 ) 


In this equation, H is the total electronic and nuclear 
Hamiltonian which is given by 


H{z,p,R;t) 


2mjv + rrie 
Am^me 
1 ^ 

mjv 


iy_ 

dp"^ p dp dz"^ 
+ Vciz,p,R;t), 


( 2 ) 


with 

Vc{z,p,R,t) = - ^ =- ^ = 

R Zuijsi T 

In these equations, Eq is the laser peak amplitude, rrie 
and rriN are masses of the electron and single nuclei, uj 
angular frequency and /(t) is the laser pulse envelope, 
which rises linearly during the first two cycles, then is 
constant for a few cycles and decreases during the last 
two cycles. Eor example, the shape of the ten-cycle laser 
pulse used in this work is shown in Eig. 1. 



FIG. 1. (Color online) The ten-cycle laser pulse shape with a 
two-cycle linear turn-on, followed by a constant intensity, and 
a two-cycle linear turn-off at 800 nm wavelength (cu = 0.057 
a.u.) and 7=4x10^^ W/cm^ intensity. 


The TDSE is solved using unitary split-operator meth¬ 
ods [30-31] with 11-point finite difference scheme through 
a general nonlinear coordinate transformation for both 
electronic and nuclear coordinates which is described in 
more details in our previous works [32-34]. The grid 
points for z, p, and R coordinates are 500, 83, and 210, 
respectively. The finest grid size values in this adaptive 
grid schemes are 0.13, 0.1, and 0.025, respectively for 2 :, 
p, and R coordinates. The grids extend up to Zmax = 98, 
pmax = 25, and Rmax = 16. The HHG spectra are calcu¬ 
lated as the square of the windowed Eourier transform of 
dipole acceleration az{t) in the electric field direction (z) 
as 

, r ,2 

S{uj) = I / < I > H{t) exp[—iujt] dt \ , (4) 

Jo 

where 

H{t) = ^[1 -cos{2tt^)], ( 5 ) 

is the Hanning filter and T is the total pulse duration 
(one optical cycle of 800 nm wavelength equals 2.6 fs). 

To obtain contributions of different vibrational and 
electronic states to total HHG spectrum, we first decom¬ 
pose the total wavefunction as [20] 

i’iz, P, R', t) = ELi Cj(-R; t)'ij:i{z, p; R) + Aesiz, p, R; t). 

(6) 
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FIG. 2. (Color online) High-order harmonic spectra produced by the NBO HJ, DJ, XJ and BO (HJ(BO)) under seven- 
cycle laser pulses of 800 nm wavelength at / =4 x 10^^ W/cm^ intensity. For better clarity, the range of 1-29 and 31-59 harmonics 
of the spectra are shown in the left and right panels, respectively. 


2 pi{z, p; Rys are four lowest Born-Oppenheimer elec¬ 
tronic wavefunctions of the system and Ci{R]t) =< 
yji{z, p; R)\yj{z, p, R;t) >z^p which is integrated over z 
and p electronic coordinates. In Eq. ( 6 ), ci(i?;t) = 
Cg{R]t), C2{R;t) = Cu{R;t), p; R)='tpg{z, p-, R) 

and p] R)='ipu{z^ P] R) in which 'ipg{z,p;R) and 

yju{z,p;R) are ground and first excited electronic wave- 
functions, respectively, corresponding to the Isag and 
2pau states. If we substitute Eq. 6 to Eq. 4 and re¬ 
tain the dominant terms, we arrive at [ 20 ] 

Stot — Sg{^) + 5'n(cj) + 2[Ag{Lu) X Au{id)], (7) 
where Sg{uj) = \Ag{uj)\‘^ and Su{oo) = \Au{uj)\‘^ and 
Ag{Uj) = 


in which S^{u) = \A^{uj)\'^ with 

Al{u) = ( 10 ) 

n—1 « 

X / 2i?e < Cu,{t)'tp„,{R)'tpg{z, p-, R) I az{t) \ tpres(t) > 

which is sum over the vibrational states with vibrational 
quantum number Ui and Cy.{t) =< 'ipy.{R)\cg{R;t) >. 
The dissociative term Sd{oo) = \Ad{oj)\‘^ with Ad{uj) = 
Ag{uj) — Ai,{uj). The last term in the Eq. (9), 2[Al{uj) x 
Ad{uj)]^ corresponds to the vibrational interference term 
(VIT) between the bound and dissociative terms of the 
Isag state. The time profile of harmonics is obtained by 
Morlet-wavelet transform of dipole acceleration az{t) via 
(8|35-36] 


j 2Re < Cg{R-,t)'tpg{z, p] R) \ a^it) \ i^res{z, p, R;t) > 
Auiuj) = 

J 2Re < Cu{R-,t)ipu{z, p; R) \ a^t) \ ipresiz, p, R;t) > 


w{uj^ t) = 


UJ 


'K‘^a 


/ +c» 

' 

-c» 


az{tyexp\—i{jj{t' — t)\exp[—- 




]dt'. ( 11 ) 


We set cr = 27r in this work as used in [36]. 


Sg{uj) and Su{oo) denote recombination to the Isag and 
2pau states, respectively and the term 2[Ag{uj) x Au{uj)] 
corresponds to the electronic interference term (EIT) be¬ 
tween these two electronic states [20]. In this work, the 
second {'ips{z^ p; R)) and third ( 7 /^ 4 (z, p; i?)) excited elec¬ 
tronic states are not considerably populated (both reach 
maximum ^ 1.4% while the population of the first ex¬ 
cited state reaches up to 40%) during the interaction and 
therefore the corresponding terms Ss{uj) and 5 ' 4 (ci;) are 
negligible and consequently have been ignored. 

We also decompose the Sg{uj) into the bound and dis¬ 
sociative parts as 

— ^b{^) + Sd{bj) + 2[Al{uj) X Ad{uj)]^ (9) 
in which the bound term Si){uj) = and Ai,{uj) = 


III. RESULTS AND DISCUSSION 

The HHG spectra of HJ, DJ, XJ and (BO) ob¬ 
tained under seven-cycle 800 nm laser pulses of 1=4x10^^ 
W/cm^ intensity are shown in Eig. 2. Eor better visu¬ 
alization, harmonic orders between 1-29 and 31-59 are 
shown separately. Eor the fixed-nuclei approximation 
case (H^(BO)), only odd harmonics are seen while for 
NBO cases, although odd harmonics are dominant at low 
orders but the HHG signals are complicated. Eor Xj 
isotopomer, small redshifts in harmonics relative to odd 
harmonics are seen which are more discernible for higher 
harmonics. Eor DJ and HJ cases, the patterns of the 
redshift are complicated and the HHG spectra broaden 
for most harmonic orders. 
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Time (laser cycles) 

FIG. 3. (Color online) The Morlet-wavelet time profiles 
for NBO H+, D+ and X+ and BO H+ (H+(BO)) under 
seven-cycle laser pulses of 800 nm wavelength and 7=4 xlO^^ 
W/cm^ intensity. The HHG intensities are depicted in color 
logarithmic scales on the right side of panels. 


To demonstrate the origin of redshifts discussed above, 
we plotted temporal dependence of harmonic orders for 
both BO and NBO cases in Fig. 3. As can be seen, 
mostly short trajectories contribute to the HHG because 
long-trajectory suppression is considerable due to the 
wavepacket spreading and also nuclear motion [21]. For 
laser pulses of sin^ and Gaussian envelopes, we can at¬ 
tribute rising and falling parts to the laser pulse. The 
blueshift and redshift of harmonics occur in rising and 
falling parts, respectively (see Ref. [12], and references 
therein). Bian et al showed that the HHG occurs mostly 
at the laser falling part because of the nuclear motion and 
enhanced ionization, giving rise to the redshift of har¬ 
monics [12]. In this work, we used trapezoidal envelopes 
having a two-cycle rising and a two-cycle falling part. As 
it is seen in Fig. 3, for HJ, DJ and Xj, there are three 
relatively intense peaks between 5 to 7 laser cycles, i.e. 
in the falling part of the laser pulse, while there is not 
any significant peak between 0-2 optical cycles for the 
laser rising part. For Xj, we see several intense peaks 


before five optical cycles and therefore the redshifts from 
laser falling part become less important in comparison 
with HJ and Dj and for Dj similarly the redshifts are 
less than those in HJ. To see whether the redshifts are 
due to laser falling part, we plotted the HHG spectra of 
three isotopomers in Fig. 4 which is similar to Fig. 2 
but here Fourier transform is integrated over the first 5 
optical cycles to remove the contribution of the last two- 
cycle falling part of the laser pulse. As it is obvious, the 
HHG spectra have been smoothed and complicated fea¬ 
tures, already seen in Fig. 2, are not seen in Fig. 4 and 
also the redshifts have been decreased. 

If we choose a longer trapezoidal laser pulse (with¬ 
out changing the optical cycles of the rising and falling 
parts of the laser pulse), we expect that contribution of 
laser falling part becomes less as the HHG process occurs 
mostly before laser falling part especially for the lighter 
isotopes with faster nuclear motions and more population 
depletion. To show this, we calculated the HHG spectra 
for ten-cycle laser pulses (as shown in Fig. 1) with other 
laser parameters similar to those of Fig. 2. The results 
are shown in Fig. 5 with their corresponding time pro¬ 
files depicted in Fig. 6. As can be seen in Fig. 5, the 
redshifts are small in comparison with Fig. 2 and for 
HJ we can say that no redshift is observed but for DJ 
and X^ small redshifts are seen. Fig. 6 shows that for 
HJ and DJ most HHG occurs before the pulse falling 
part. The intensity of harmonic peaks decreases for H^ 
and DJ at time t > 6.5 o.c. and at t > 7.5 o.c., respec¬ 
tively. The population of H^ (Dj) is depleted from 0.8 
at t = 6.5 o.c. (7.5 o.c.) to 0.15 (0.45) at the end of the 
ten-cycle laser pulse . Therefore, these population deple¬ 
tions are responsible for the suppression of the HHG at 
the last two-cycle laser pulse in Fig. 6. For X^, the pop¬ 
ulation depletion is minimum due to slow nuclear motion 
not allowing the nuclei reach large internuclear distances 
where enhanced ionization occurs. Therefore the HHG 
has considerable intensity in the falling part of the laser 
pulse. Fig. 6 shows that for Xj there are three HHG 
peaks having redshift at the falling part of the laser pulse 
(between 8-10 optical cycles) and there are twelve HHG 
peaks without redshift before the falling part (between 
2-8 optical cycles) which reduce and cover the contribu¬ 
tion of the three peaks with redshift into the total HHG 
spectrum as we do not observe any redshift in the total 
HHG spectrum of Xj in Fig. 5. In Fig. 7, the ion¬ 
ized outgoing wavepacket (10W) probability density in 
the 2; direction as a function of internuclear distance for 
the ten-cycle laser pulse is shown for H^, and X^. 
This figure shows that the ionization occurs for DJ and 
HJ mostly at large internuclear distances where charge- 
resonance enhanced ionization happens. It must be noted 
for Fig. 7 that we calculated the lOW probability density 
from the boundary on the electronic z coordinate by the 
virtual detector methods [32]. The nuclei in this lOW 
are mainly under influence of the Coulomb repulsion of 
each other. Therefore, we expect the internuclear dis¬ 
tance of the lOW increases between the ionization time 
and the time to reach the virtual detector. This increase 
of internuclear distance is more for a smaller internuclear 
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FIG. 4. (Color online) The same as FIG. 2, but Fourier transform is done over first five optical cycles (set T = 5 o.c. in Eq. 

(4)). 


CO 

o 

D) 

O 



Harmonic Order Harmonic Order 


FIG. 5. (Color online) The same as FIG. 2, but for a ten-cycle laser pulse shown in Fig. 1. 


distance. Therefore, if we want to introduce the lOW 
probability density in Fig. 7 as the ionization probability, 
it is relatively deformed and shifted to larger internuclear 
distances. 

The time-dependent population as function of internu¬ 
clear distance is also shown in Fig. 8 for the H^, 
and Xj for the ten-cycle laser pulse. The probability of 
finding the molecule at larger internuclear separations in¬ 
creases in time which is more for the lighter isotopomers. 
The intensity decrease of harmonic peaks in Fig. 6 for 
H+ {t > 6.5 o.c.) and Dj (t > 7.5 o.c.) is consistent with 
finding the molecule at large internuclear distances (Fig. 
8) where ionization (Fig. 7) and therefore the population 
depletion is maximum. 

In order to figure out the contributions of different elec¬ 
tronic and vibrational states into the redshifts observed 
in Fig. 2 for the seven-cycle laser pulse, we plotted sepa¬ 
rately recombination to dominant electronic states, 
and 2p(Ju^ and different vibrational states (V^zyJ of 
state which are depicted in Figs. 9-11 for X^, and 
HJ, respectively. We should note that for Figs. 9-11 
we have depicted harmonic orders greater than 19 to fo¬ 
cus mainly on the recombination of continuum electron 


into bound states. Transitions between bound electronic 
states also lead to generation of harmonic orders up to 
ionization potential Ip (in this work Ip = 1.1 correspond¬ 
ing to harmonic order ^ 19), leading to further complex¬ 
ity of harmonic orders below Ip which was not intended 
in this work. In Fig. 9(a) for Xj, the spectra of Sg + Su 
and Slot have nearly the same magnitude and thus we can 
conclude that EIT has not significant effect. It is also ob¬ 
served for recombinations to Iscr^ and 2pau states that Sg 
is dominant over Su for almost high harmonic orders but 
at low harmonic orders Su is slightly higher (Fig. 9(b)). 
As Sg is dominant for most harmonic orders, so we focus 
on its components, 5'^ and Sd^ to find out their contribu¬ 
tion in the redshifts seen in Fig. 2. As can be seen in Fig. 
9(c), Sd is much lower than S^ and consequently disso¬ 
ciative part of Sg has insignificant portion. Comparison 
of Sg and Sh + Sd in Fig. 9(b) indicates that the VIT is 
negligible. Therefore vibrational components of 5'^ should 
be the source of redshifts observed for Xj. To find out 
how many vibrational states contribute to 5'^, we com¬ 
pared different 5'^ with Sh and found that Sl^ matches 
well with the S'^, i.e. mainly up to 10 vibrational states 
contribute to the Sb- Therefore, lowest vibrational states 
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FIG. 6. (Color online) The same as FIG. 3, but for a ten-cycle 
laser pulse. 


{^i <9) mainly contribute to the Stot and the higher-lying 
vibrational states {ui >9) have little contribution. It was 
expected because of the heavy nuclei of X J which makes 
nuclear motion slower and thus the nuclear wavepacket 
can not reach large internuclear distances at these seven- 
and ten-cycle laser pulses and therefore high vibrational 
states are not significantly populated. We also already 
observed that Sd (Fig. 9(c)) was also very weak which is 
consistent with contribution of the lowest bound vibra¬ 
tional states discussed above. If we compare and 
in Fig. 9(d), the redshifts are larger for the demon¬ 
strating that the redshifts in the HHG spectra originate 
from the contribution of higher vibrational states of Isag 
state. The patterns of the HHG spectra for to Sl^ 
is similar to those of and only the HHG signals have 
been intensified, with no extra shift compared to what 
observed for S^. 

Now we investigate contributions of decomposed com¬ 
ponents of total HHG spectrum for in Fig. 10 as 
applied for the Xj. The DJ is much lighter than Xj 
and therefore we expect a faster nuclear dynamics and 
different coupled electron-nuclear dynamics. By compar¬ 



Internuclear distance 


FIG. 7. (Golor online) The ionized outgoing wavepacket 
(lOW) probability density in the z direction as a function of 
internuclear distance for NBO HJ, DJ and Xj under ten-cycle 
laser pulses of 800 nm wavelength and 7=4 xlO^^ W/cm^ 
intensity. 


ing Sg + Su and Stot ia Fig. 10(a), it can be deduced 
that the FIT is not important for low harmonic orders 
but has a pronounced effect on high ones (>39), lead¬ 
ing to considerable modulations on the signal. The EIT 
weakened the signal, partially around odd-order harmon¬ 
ics and remarkably around even-order harmonics. In Fig. 
10(b), for low harmonic orders (<39), the intensities of 
Sg and Su differ considerably which is the reason why 
the EIT is insignificant at these harmonic orders. But 
for harmonics greater than 39, Sg and Su have almost 
close signal magnitude and this is while, as mentioned 
above, the EIT is more noticeable for these harmonic or¬ 
ders. One of the important points observed for Sg and 
Su is that the redshifts occur for both components. The 
2pau state is an anti-bonding (dissociative) state and this 
pure dissociative electronic state has also contributed to 
the observed redshifts. Similar to X^, Sd is much less 
intense than Sb and thus its contribution to Sg is negli¬ 
gible. Comparing Sg with Sb + Sd in Eig. 10(b) reveals 
that the VIT is insignificant similar to what observed for 
X 2 • Eig. 10(d) shows that the number of vibrational 
states contributing to the Sb has been increased in com¬ 
parison with Xj. The bound vibrational stationary states 
of the electronic ground state of Xj are more and more 
compact than those of DJ, thus contribution of more vi¬ 
brational states of the latter into the HHG indicates that 
the DJ has reached larger internuclear distances. That 
means that vibrational states with higher energies (wider 
nuclear distribution) are populated for the D J and there¬ 
fore contribute to the HHG. Eig. 10(d) shows that with 
increasing the number of z^^’s in Sb (Eq. 10), the redshift 
of odd harmonics has become larger and even harmonic 
orders has gotten more intense due to the contribution 
of higher vibrational states. Contribution of high vibra¬ 
tional states to the formation of even harmonic orders has 
been also recently reported for a one-dimensional H^ in 
relatively long laser pulses (10- and 14-cycle pulses) [18]. 
We mentioned above that the EIT has dominant effect for 
high harmonic orders (Eig. 10(a)). Another point that 
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FIG. 8. (Color online) Time-dependent population for NBO 
HJ, DJ and Xj as function of internuclear distance under 
ten-cycle laser pulses of 800 nm wavelength and /=4 xlO^^ 
W/cm^ intensity. 



Harmonic order 

FIG. 9. (Color online) The components of the HHG spec¬ 
trum of the X .2 exposed to seven-cycle laser pulses of 800 nm 
wavelength and /=4 xlO^^ W/cm^ intensity. 


we can add about the EIT is that its effect, as already 
mentioned, is seen to be more noticeable around even 
harmonic orders rather than those close to odd harmonic 
orders. As contribution to even-generated harmonic or¬ 
ders comes mainly from high-lying vibrational states and 
these states are closer in energy to the 2pau state than low 
vibrational states, the EIT affects considerably the HHG 
signal around even harmonic orders than those close to 
odd harmonic orders. 

Most above justifications for observations in D J can be 
considered for HJ and we just note some specific features 
of the latter. The EIT has increased the signal for low 
harmonic orders but has decreased it for high ones (Eig. 
11(a)). The signal Su is higher than Sg for harmonic or¬ 
ders 39-53 and the redshifts are seen for both Sg and Su 
in Eig. 11(b). Due to faster nuclear motion of H^ rela¬ 
tive to DJ and Xj, the molecular ion can reach large in¬ 
ternuclear distances in which charge-resonance enhanced 
ionization occurs and leads to considerable population in 
2pau state. Unlike X^ and D^, Sd shows comparable 
intensity as Sb and even for some high harmonic orders 
is higher than Sb as shown in Eig. 11(c). Comparing 
Sg with Sb + Sd in Eig. 11(b) shows that the VIT could 


change (weaken or enhance) the HHG spectrum not ob¬ 
served for DJ and X^. In Eig. 11(d), as the number of vi¬ 
brational states that contribute to 5'^ increases, the peaks 
close to odd harmonics are slightly being red-shifted and 
the HHG spectrum near even harmonics get more intense. 

Eor the case of the ten-cycle laser pulse, all the HHG 
components, Stot^ Sg^ Su, Sb, Sd and 5'^’s show only 
odd harmonic orders with very small redshift compared 
to those observed for the seven-cycle laser pulse in Eig. 
2. The contribution of the above-mentioned decomposed 
components for each isotopomer is similar (only around 
odd harmonic orders) to those observed and mentioned 
in Eigs. 9-11. The only difference seen is that the VIT 
is insignificant for the ten-cycle laser pulse case, even for 
the H+. 


IV. CONCLUSION 

We solved numerically full-dimensional electronic time- 
dependent Schrodinger equation for HJ and its iso- 
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Harmonic order 

FIG. 10. (Color online) The same as FIG. 9, but for DJ. 


topomers with and without Born-Oppenheimer approxi¬ 
mation under seven- and ten-cycle trapezoidal laser pulses 
of 800 nm wavelength and 7=4x10^^ W/cm^ intensity. 
The effects of the falling part of trapezoidal laser pulses 
and control of these effects on the HHG spectrum were 
investigated by considering different pulse durations and 
isotopomers. 

In the case of the fixed-nuclei calculation, only odd har¬ 
monic orders are generated but for the case of the freed 
nuclei (beyond the Born-Oppenheimer approximation), 
the HHG spectra change depending on the isotopomer. 
The results show that considerable redshifts appear for a 
shorter laser pulse (seven cycles) in the presence of nu¬ 
clear motion, originating mainly from the last two-cycle 
falling part of the laser pulse. Whereas, by considering 
a longer laser pulse (ten cycles), the redshifts and com¬ 
plexity of the HHG spectra are considerably reduced and 
controlled because the HHG at the falling part of the laser 
pulse has a little contribution on the total HHG spectrum. 
For DJ and HJ under ten-cycle laser pulses, the HHG is 
suppressed considerably at the laser falling part due to 
the population depletion. 

In order to get insight into the underlying components 
of the complicated patterns observed for the seven-cycle 
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Harmonic order 


FIG. 11. (Golor online) The same as FIG. 9, but for HJ. 


trapezoidal laser pulses, the HHG spectrum was decom¬ 
posed into different components representing recombina¬ 
tions to different electronic states and 2pcr^. We 

also decompose recombination to ground electronic state 
^s(jg into the bound {Sh) and dissociative {Sd) compo¬ 
nents. For the heaviest isotopomer, Xj, mostly ground 
electronic state and low vibrational states contribute to 
the HHG spectrum and high harmonic orders show a red- 
shift coming mainly from the high vibrational states of 
the ^-SCFg state. For isotopomers Dj and Hj, recombi¬ 
nation to the first excited electronic state 2p(Ju becomes 
also significant and the interference term (EIT) between 
^s(jg and 2p(Ju states plays an important role and changes 
the HHG signal, remarkably around even-order harmon¬ 
ics. More vibrational states contribute to the HHG signal 
for lighter isotopomers indicating that larger internuclear 
distances have become accessible during the interaction 
due to the faster nuclear motion. With increasing the 
number of vibrational states, the HHG signal of odd har¬ 
monic orders is being slightly redshifted and that of even- 
order harmonics gets more intense. The dissociation part 
Sd and VIT seen to be noticeable for the HJ. 

Another way that we can suppress the effect of the 
falling laser pulse on the HHG spectrum is to use enough 
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high laser peak intensities, instead of a longer laser pulse, which should be considered, 
not considered in this work. We should note that, at 
very high laser intensities or longer laser rising parts, the 

blueshift of harmonics at the laser rising part may arise V. REFERENCES 
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